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We present a detailed account of the two-photon states generated by SPDC in both type I and type 
II phase matching, including the effects of anisotropy of the nonlinear medium and the frequency 
spread of the down-converted fields. Accurate as well as simplified expressions are derived for type 
I and type II phase matching in the context of Fourier Optics. The main results are compared with 
experimental data available in the literature, showing good agreement in all cases. 
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I. INTRODUCTION 

In the course of the last three decades, spontaneous 
parametric down-conversion (SPDC) has proven to be a 
valuable tool in the experimental investigation of fun- 
damental properties of the electromagnetic field in the 
quantum domain [l|, 0, 9 , including nonclassical correla- 
tions, entanglement and nonlocality. The fact that SPDC 
is capable of generating pairs of photons in a wide range 
of frequencies and wave vectors, in addition to the fact 
that these photons may be entangled in a number of dif- 
ferent degrees of freedom, qualifies SPDC also as a unique 
tool in the demonstration of quantum information pro- 
cedures and protocols [1]. 

In many applications, especially in the pioneering ones, 
a basic knowledge of a few properties of the two-photon 
states generated by SPDC is enough to explain the ef- 
fects discussed. As the applications become more sophis- 
ticated, for example in those that combine correlations in 
more than one degree of freedom, a more comprehensive 
approach is necessary in order to enhance the capabilities 
of quantum state engineering. In this direction, the use of 
Fourier Optics concepts in two-photon optics represents 
a significant step. 

An important application field of SPDC is the con- 
struction of entangledphoton state sources with high 
fidelity and purity @, i, 0, H, !, M, El El- In order 
to obtain efficient and dependable two-photon sources, a 
number of questions must be addressed, as for example, 
phase compensation, mode coupling, pair collection effi- 
ciency, maximum photon flux, etc. To be optimized, all 
these points require a detailed knowledge of the SPDC 
process. 

Another context in which a more precise model is cru- 
cial refers to the issue of conservation of orbital angu- 
lar momentum (QAM) of light in the process of SPDC 
0, El, El El El 111- Oversimplifications or incon- 
sistent assumptions about the two-photon state gener- 



ated have led to conflicting conclusions on whether OAM 
is conserved or not in SPDC. 



In this work, we present a detailed account of the two- 
photon states generated by SPDC in both type I and type 
II phase matching, including the effects of anisotropy of 
the nonlinear medium and the frequency spread of the 
down-converted fields. As an exact approach to this 
problem presents a considerable level of difficulty, we are 
forced to restrict ourselves to some approximations. We 
adhere to the commonly adopted simplified procedure for 
the quantization of the electromagnetic field in the non- 
linear medium by just multiplying the k vectors by the 
corresponding refractive indices and keeping the vacuum 
expressions for the field operators. The frequency spread 
of the down-converted fields is considered to be small 
compared to their central frequencies . We also work 
in the paraxial approximation, in which the propagation 
of the fields is much easier to deal with. Fortunately, 
these approximations, are good enough to encompass a 
great number of practical cases. We admit, however, that 
while the model discussed here is helpful in far field appli- 
cations, it may not be appropriate in discussions of the 
fundamental interactions inside the nonlinear medium. 
Whenever possible, the results reported here are com- 
pared with experimental data, obtained in our labs or 
reported in the literature. 



It should be mentioned that other authors have re- 
ported their contributions to this field, in different levels 
of detail and scope [H, [H, HH, H3, [H, [H, HE [H, H3, 
|28| . From our point of view, a comprehensive descrip- 
tion of SPDC that accounts for all features observed so 
far in two-photon states, with expressions ready to use in 
Fourier optics, is still missing. It is in this context that 
we present our work. 
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II. TWO-PHOTON STATE GENERATED BY 
SPDC 



The interaction Hamiltonian describing the optical 
processes of parametric down-conversion in nonlinear 
birefringent crystals is usually written in terms of a sim- 
plified field quantization in matter [291 ] . The field modes 
participating in this process are coupled by the second- 
order susceptibility tensor ■ Due to the phase match- 
ing conditions, the coupling occurs only for some sets of 
polarizations, known as type I and type II. In negative 
uniaxial crystals (n e < n ) such as BBO (/?— BaB 2 C>4) 
and Lithium Iodate (L1IO3), the down-conversion pro- 
cess with type I phase matching can be summarized as 
e — * 00, meaning that one photon of the pump beam with 
extraordinary polarization is converted into two photons 
with ordinary polarization. In type II phase matching, 
the process is represented by e — > oe or e — > eo. 

In a perturbative approach, the (post-selected) two- 
photon state generated by SPDC is written as [2(| 

l*> = X) H $ k pkl k 2CTpCTlCT2 |ki,(ii) |k 2 ,a 2 ) , (1) 

01,02 ki ,k 2 

where ki and k 2 are the wave vectors of the down- 
converted fields, k p is the wave vector of the pump field, 
whose frequency lu p is assumed to be well defined. <jj 
indicates the polarization of each down-converted field, 
that can be ordinary (0) or extraordinary (e), and |kj, <jj) 
represent a one-photon state in the plane wave mode 
kj , <jj . Considering that the nonlinear crystal is a rect- 
angular block of sides L x , L y , L z , with two faces parallel 
to the plane xy, the amplitude $ is given by 



and L02 can therefore be written as 



{k,<r} 



t e { 2> sine 
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x £k p ,<r p / ' 



dr, 



(2) 



where {k, a} represents the set of indices kpkik 2 crp(7i<7 2 , 



, V is the quantization volume, 



n(\ij, Uj) is the 



ffk 3 -,u 3 = " y2e Vn 2 (\t ] ,<7 j ) 

X is the interaction volume (L x L y L 
refractive index corresponding to the mode kj ,aj, f2 = 

(2) 

u>i + L02 — u) p , t is the interaction time, xljk IS related 
to the second-order nonlinear susceptibility tensor [2^ |. 
(eki,<Ti)j (j — x, y, z) are the cartesian components of the 
polarization vectors, £k p ,<r p is the pump field amplitude 
in the mode k p , a p , and A = ki + k 2 — k p . 

In order to simplify expression it is convenient to 
make the following approximations: 

(a) The pump frequency uj p is well defined and the 
interaction time is long, so that the term sincf2T/2 is 
significant only when uj\ + c<j 2 = lu p . The frequencies u>\ 



(1 



-2 = -yd 



(3a) 
(3b) 



This assumption can be justified by the use of a moder- 
ate power continuous-wave pump laser so that the time 
interval between two down-conversions is large compared 
to the detection resolving time. 

(6) The frequency spread of the detectable down- 
converted fields is small compared to the central fre- 
quency ( \v\ -c 1), so that the dispersion of the refractive 
indices around the central frequency uj p /2 is small and 
a linear approximation can be used. This assumption is 
justified by the use of narrow-band interference filters in 
front of the detectors. For example, for a lOOnm wave- 
length spread centered at 700nm, v lies in the interval 
(-0.08,0.08). 

/ 2 \ 

(c) The terms g^^j and x\jk are slowly- varying func- 
tions of kj, so that they may be taken as constants in 
the intervals considered for kj . 

(d) The pump beam propagates along the z axis and 
the crystal is large enough in the x and y directions to 
contain the whole pump beam transverse profile. In this 
case, L x and L y can be extended to infinity and the last 
term in expression ([5]), the integral JV e is pro- 
portional to 



5(k lx + k 2 ) 5(kiy + k>2y - k pv ) 



z c +L z /2 



e~ lA * z dz, 



-L z /2 



where z c locates the center of the crystal. 

(e) The quantization volume is large enough to justify 
the replacement of summations in k by integrals. 

(/) The pump beam contains only extraordinary po- 
larization. It is implicit in this assumption that we are 
dealing with negative birefringent crystals. 

Under the above assumptions, Eq. ([1]) is written as 



I*) = j dl, J dc t y~ d q2 < i , ( T lCT2 (qi,q2,^) 



£71 -CT2 

x \qx,i/,ax) |q 2 , ~v,a 2 ) 



(4) 



where \(\j,v,Uj) represents a one-photon state in the 
mode defined by the transverse (xy) component of 
the wave vector, the frequency (1 + v)to p /2 and by the 
polarization Uj . The amplitude $ is now reduced to 

$ CT1CT2 w C aia2 G(v) £ (qi+q 2 ) / C+ *' e- lA ' z dz, (5) 

Jz c -Lz/2 

where C aia2 is a coupling constant, which depends on the 
nonlinear susceptibility tensor, G(v) is the spectral func- 
tion defined by the narrow bandwidth filters placed in 
front of the detectors, and £ is the plane wave spectrum 
of the pump beam £(q p ), with q p replaced by qi + q 2 . 



The integral in z is discussed in what follows. The func- 
tion G(y) may be centered at zero in frequency degener- 
ate configurations or at some fixed small detuning vq. 

Most of the properties of the down-converted fields is 
determined by the longitudinal wave vector mismatch 



A 2 = ku + k 2z - k 



)>--• 



(0) 



If the anisotropy of the medium is neglected (which is not 
always convenient), v = (monochromatic approxima- 
tion), and the crystal is cut for collinear phase matching, 
A z has a simple expression in terms of qi and q2, with 
the z components replaced by kj z = \/\kj\ 2 — ]0i] 2 ' In 
the paraxial approximation, kj Z can be approximated by 
k jz w |kj| - Iq^ | 2 /2|k J |. This leads, for z c = 0, to 



+ W2 / I 

e -tA„z ^ z K ginc 



4|k K 



|qi - 02 1 



which appears in the amplitude reported in Ref. [30J. 
This result, however, is a good approximation only when 
the crystal is very thin, the diffraction (Rayleigh) length 
z of the pump beam is large, and the detection area is 
small so that the sine function can be approximated by f . 
In more general conditions, the anisotropy of the crystal 
must be taken into account. 

From this point on, we will refer to the arrangement 
depicted in Fig[TJ 




FIG. 1: Geometry of the problem. The uniaxial medium is 
assumed to be a slab of thickness L z , having its faces parallel 
to the xy plane. The optic axis lies in the xz plane, at an 
angle 6 with the z direction. 



7 = 



n 2 sin 2 9 + n 2 cos 2 I 



(11) 



(12) 



with 



III. THE EFFECT OF ANISOTROPY 

Let us consider a monochromatic electromagnetic held 
in a plane wave mode e I ( k - r ~ wt ) propagating through 
a nonmagnetic uniaxially birefringent medium with or- 
dinary and extraordinary refractive indices n Q and n e , 
respectively. Considering that the optic axis lies on 
the plane xz, making an angle 9 with the z direction, 
and that the plane wave has extraordinary polariza- 
tion, Maxwell's equations require that the components 
of the wave vector k, whose cartesian components are 
(q x ,q y ,k z ), satisfies the ray surface equation [3l[ 



(n 2 cos 2 9 + n 2 sin 2 9)q 2 x + (n 2 a sin 2 9 + n\ cos 2 9)k 2 z + 
2(n 2 -n 2 e ) sin 9 cos 9 q x k z +n 2 q 2 y = [n n e ^j . (7) 

Solving Eq. for k z , we find 



where 



k z = -aq x + J k 2 - Peg. - -yq 2 ,, 



{n 2 — n 2 ) sin 9 cos 9 



i 2 , sin 2 9 + nl cos 2 9 ' 



= 



n n e 



n 2 , sin 2 9 + n 2 cos 2 9 J 



(8) 



(9) 



(10) 



rj = 



n 2 sin 2 9 + n 2 cos 2 9 



(13) 



The ordinary and extraordinary refractive indices n and 
n e are obtained from the Sellmeier equations. For BBO, 
n and n e are obtained from [32J 



n% = 2.7359 
n 2 = 2.3753- 



0.01878 
A 2 -0.01822 
0.01224 



0.01354 A 2 , 
0.01516 A 2 , 



A 2 -0.01667 
with A given in /zm. For Lithium Iodate we use [3. 

1.332068 A 2 



2.083648- 



1.673463 



A 2 - 0.035306 
1.245229 A 2 



A 2 - 0.028224 
Within the paraxial approximation, 

1 



0.008525 X, 
0.003641 A 2 . 



aqx 



2k 



{0q£ + 7<z 



(14a) 
(14b) 

(15a) 
(15b) 

(16) 



The term a represents a linear displacement in the x 
direction, also known as the walk- off. It reaches its max- 
imum value in the neighborhood of 9 = 45°. The terms (3 
and 7 account for the deviation of the curvatures of the 
ray surface from a spherical surface in the x and y direc- 
tions, respectively. Differently from a, the terms (3 and 7 
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FIG. 2: (Color online) Values of a for for 300nm and 600nm in 
BBO and Lithium Iodate as functions of the phase matching 
angle 9. The curves are ordered from top to bottom, according 
to the legend. 



FIG. 4: (Color online) Values of 7 for for 300nm and 600nm in 
BBO and Lithium Iodate as functions of the phase matching 
angle 9. The curves are ordered from top to bottom, according 
to the legend. 
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FIG. 3: (Color online) Values of for for 300nm and 600nm in 
BBO and Lithium Iodate as functions of the phase matching 
angle 9. The curves are ordered from top to bottom, according 
to the legend. 



have little influence in the phase matching. The term 77 is 
the refractive index for a plane wave with extraordinary 
polarization, propagating along the z direction, whose 
wavenumber is x. Figs. [2j [3] and [4] show the values of a, 
(3 and 7 as functions of the phase matching angle 9 for 
300nm and 600nm in BBO and Lithium Iodate. 



IV. TYPE I PHASE MATCHING 

Let us consider first the case of type I phase match- 
ing, where the pump field has extraordinary polarization 
and the down-converted fields have ordinary polarization 
(e — * 00). In this case, 



1 

2 k, 



2 

px 



■ Iplly) 



(17) 



hz ~ fa - 



2V 



«22 



kn- 



2k 2 



(18) 



(19) 



where fa — n i u>\/c and fa — n o2 ll>2 / c. Since we are 
assuming v -C 1, the refractive indices n \ and n 2 can 
be written in a linear approximation as 



n 2 



n (l + av), 
n (l - av), 



(20a) 
(20b) 



where n a is the ordinary refractive index at the frequency 
u>p/2, and factor a is given by 



uj p dn 
2n n doj 



(21) 



j=uj p /2 



In order to calculate A z , which in type I phase match- 
ing will be referred to as A°°, we have to write k\ z + 
k2z - k pz replacing q px and q py in Eq. (JTTJ) by q lx + q 2x 
and qi y + q 2yi respectively. After some manipulation, we 
find 

A°° « K + a p (q lx + q 2x ) 
1 



n D K 
b 



^{qix + Q2x) 2 



9 1 

2 (?ly 



(22) 



where K = lu p /c, b = l3 p n /i] p , and g = J p n /r)p. The 
subscript p means that a, (3 and 7 refer to the pump 
field. Hoo is the type I collinear index mismatch, defined 
as 



fJ-o 



n (l 



av 



(23) 
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FIG. 5: Angular components of the down-converted fields. 



In the derivation of Eq. (|22|) . terms in v 2 were ne- 
glected with respect to v. For type I collinear (/i QO = 0) 
phase matching in BBO pumped with X p — 351nm, 
a p = 0.0747, a = 0.02, b = 1.06, and g = 1.11. For 
Lithium Iodate, a p = 0.0871, a = 0.03, b = 0.951 and 
g = 1.07. 

Sometimes, instead of the transverse components of 
the wave vectors, one is interested in the output angles 
of the down-converted photons. In general, the output 
fields of interest propagate at small angles with respect 
to the z axis, so that one can approximate the output 
angles by 



c_ _ 2qxj 

Kl ' J ~ UJ^ ~ (1 + !/)«" 
t C 2 l2j 

UJ 2 (1 - v)K 



(24a) 
(24b) 



where j = x, y, and we define the vectors £ ± = (£i x ,£iy) 
and £ 2 = (£,2x,£,2y) as shown in Fig. [5] In the far field, 
£ 1 and £ 2 can be written in terms of the transverse co- 
ordinates p 1 and p 2 of the detectors as 



Ci » — Pl , 

ZD 



(25a) 
(25b) 



where zd is the z coordinate of the detection plane (par- 
allel to the xy plane). 

In terms of the angular variables £j and £ 2 > the ampli- 
tude <£> 00 , defined in Eq. (JSJ), is written as 



NG{v)E 
~KL 



y(£i + £ 2 ) + ^(£i-£ 2 ) 



xsmc 



where N is a normalization constant and 



/o (£i,£ 2 ,^) = Moo + y [(i + + (i - 



J_{(i + „)£ + (!-„)£ 



4n G 

^[(l + ^ix + (l-^)6x] 2 
2 



[(l + i ,)£ l3/ + (l (27) 



The single count rate as a function of £ 1 (£ 2 ) and i/ 
can be obtained from Eq. (f2"6")) by integration in £ 2 (^jj. 
To show that Eqs. (|26|) and (|27|) describe the correct 
dependence of $ on frequency and output angles, the 
output angle £ x is plotted as a function of the down- 
converted wavelength in Fig. [6] for a 15mm- long Lithium 
Iodate crystal pumped by a 325nm laser beam for three 
different values of the phase matching angle 9. Results 
for 9 = 59.217° (collinear) and 9 = 59.185° are in good 
agreement with experimental data reported by Bogdanov 
et al. 01 for 9 = 59.22° and 9 = 58.97°, respectively. 
The difference of 0.22° in the noncollinear phase match- 
ing angle is probably due to the fact that the noncollinear 
condition is obtained by tilting the crystal. 

The two-photon frequency spectral content of $ OD for 
a given pair of selected output angles £° and £ 2 can also 
be derived from Eq. (f2"6")l . Experimentally, this spec- 
trum can be directly measured [251 ] . or inferred from the 
Hong-Ou-Mandel dip when the measurement is taken 
with broadband filters. The expected profile of the co- 
incidence detection probability as a function of the path 
length difference 51 = cSt in a Hong-Ou-Mandel interfer- 
ometer is given by 35] 



P c (St) 



J dr 4>* (t) </>(t — St) 

1 JdTWF) 



(28) 



where, for a pump beam with a symmetric transverse 
profile, 



r)= / dv$oo&,&v) 



-iKvcSr 



(29) 



x exp [-iKzc / o0 (£i,£ 2 i v )] 



(26) 



The integrals are taken from — oo to oo. Fig. [7] shows a 
plot of P c {8t) calculated with Eqs. (J2U) and |2SJ, for a 
BBO crystal with L z = 1mm in type I phase matching, 
pumped by a 351nm laser, without interference filters. 
The collection angles are such that both down-converted 
beams are centered at 702nm. 

In the collinear detection configuration (^ 1 = £ 2 = 0), 
the spectral profile of the two-photon state depends on 
the phase matching angle 9. Figures [S] and [5] show the 
down-converted wavelengths for collinear type I phase 
matching as a function of 89, the deviation from 9 m = 
51.704° for Lithium Iodate and from 9 m = 33.543° for 
BBO, with L z = 5mm and L z = 2mm, both crystals 
pumped with A p = 351nm. 
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FIG. 6: Output angle in the xz plane as a function of the 
down-converted wavelength for a 15mm-long Lithium Iodate 
crystal pumped by a 325nm laser beam, for different values 
of the phase matching angle 8. These results were obtained 
from Eqs. (|26[) and (|27[) integrating in £ 2 with £i s = 0. The 
(nonlinear) horizontal scale was adjusted to allow a direct 
comparison with experimental data reported by Bogdanov et 
al. 13411. 
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FIG. 7: Hong-Ou-Mandel dip for a lmm-long BBO crystal cut 
for type I phase matching, with 8 = 34° pumped by a 351nm 
laser, calculated with Eqs. (|26[l and ((28}. It is assumed that 
no interference filters are used. 
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FIG. 8: (Color online) Down-converted wavelengths for 
collinear (£ 1 — £ 2 — 0) type I phase matching as a function 
of 58, the deviation from 9 m = 51.704° for Lithium Iodate 
with L z — 5mm (left) and L z = 2mm (right), pumped with 



floo = 0. 



351nm. The solid line corresponds to the condition 



L z = 5mm 



2mm 



CD 
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FIG. 9: (Color online) Down-converted wavelengths for 
collinear = £ 2 = 0) type I phase matching as a function of 
58, the deviation from 9 m — 33.543° for BBO with L z = 5mm 
(left) and L z = 2mm (right), pumped with A p = 351nm. The 
solid line corresponds to the condition jj, oo —0. 



Another useful representation of $ 00 is obtained when 
the following new coordinates are used: 



£s = 2^1+^). 

%d = 2 (^i — Ca)- 
Now, Eq. (Po| is written as 

$ 00 w NG(i>)S(Kt a + Kvt d ) 
xsinc F 00 {£ s ,£ d ,v) 
x exp [~iKz c F 00 (£ s ,£ d ,v)] , 
where N is a normalization constant and 



(30a) 
(30b) 



-2^[(l-6)€« + (1-5)4 
+(2 - b)v£sMx + (2 - g)v£ S y 



(31) 



(32) 



7 



where v 2 was neglected with respect to v. 

Many of the interesting features of $ OQ are retained 
when the anisotropy is neglected in second order, that is, 
when the parameters a and b are both approximated by 
1. In this case, 



1 



/'o 



In the frequency degenerate case (y = 0) 
Foo(€ s , £<j, v) w (i 00 + a p £ sx - 



^d 

2n n 



(33) 



(34) 



Let us analyze the behavior of Eqs. ([311 and ([34| in 
two different measurement schemes: (a) When the detec- 
tors are scanned in opposite directions in the far field, so 
that £ s = 0. The probability of coincidence detection is 
proportional to 



\®oo\ 2 d 



NG{v)£(0) sine 



KL Z 
4n n 



(R 2 - &) 



(35) 



Eq. (J35 
with a radius 



describes a circular profile in the variable £d, 
= ^/'ln [ioo and a half- width 5R = 
y/ R 2 + Aitriol KL Z — R. Under these detection condi- 
tions, the coincidence transverse profile maps the well- 
known down-converted light cones. 

(b) When the detectors are scanned in the same direc- 
tion in the far field, so that = R. This condition 
means 



(36) 



Eq. (|36p illustrates the possibility of transferring the an- 
gular spectrum £ (q) from the pump beam to the down- 
converted field £(K£ S ) = £(qx + q 2 ) [HJ, and its depen- 
dence on the crystal length L z . The sine function has 
a half- width <5£ sx = 2-k j KL z ol p in the x direction. De- 
pending on the spatial bandwidth of the pump beam, this 
factor may limit the transfer of its angular spectrum to 
the two-photon field. Considering that the pump beam 
is Gaussian, the transfer will be satisfactory only when 
its waist wq is much larger than L z a p /2i:. For a BBO 
of L z ~ 1mm, L z a p /2n ~ 10/xm. This "clipping" in 
the angular spectrum transfer, caused by the walk-off 
term a p ^ sx , has consequences in the anisotropy of en- 
tanglement in spatial variables of the two-photon state 
(36| . Figs. [TOl and fTTI show a comparison between ex- 
perimental results obtained in our labs and the predic- 
tions of Eq. (|36p for a BBO crystal with L z = 5mm 
in collinear type I phase matching, pumped by a 405nm 
laser with wo w 25^m. In Fig. [101 the two detectors are 
scanned in the same sense, along the x direction. The 
effect of the walk-off term is evident. The dashed line 
shows the expected angular profile of the pump beam. 
In Fig. [IT] the two detectors are scanned in the same 
sense, along the y direction. In both cases, interference 
filters of FWHM=10nm were used in front of the detec- 
tors. 
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FIG. 10: (Color online) Comparison between experimental 
results (•) and the prediction of Eq. (|36|) (solid line) for a 
BBO crystal with L z — 5mm in collinear type I phase match- 
ing, pumped by a 405nm laser with wo ~ 25/im. The two 
detectors are scanned in the same sense, along the x direc- 
tion. The dashed line shows a gaussian profile corresponding 
to the pump laser beam. 
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FIG. 11: (Color online) Comparison between experimental 
results (•) and the prediction of Eq. (|36p (solid line) for a 
BBO crystal with L z = 5mm in collinear type I phase match- 
ing, pumped by a 405nm laser with wo ~ 25/xm. The two 
detectors are scanned in the same sense, along the y direc- 
tion. The Gaussian predicted by Eq. (|36p coincides with the 
expected profile of the pump beam. 



V. TYPE II PHASE MATCHING 

In type II phase matching, the calculation is similar, 
except for the fact that one of the down-converted beams 
has extraordinary polarization (e — > oe or e — > eo), that 
is, 

I*) = Jdv j "(fqi Jdq 2 [® oe |qi, v\ o) |q 2 , -v; e) + 

$ eo \qi,v,e) |q 2 ,-^;o)]. (37) 
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Now, the k z component of the e-polarized down- 
converted wave has to be written according to expression 
(fl"6|) . that is, for the case e — ■> oe, 



^pz 



oQpx „ (ftpQpx + lp1py)i 



kiz ~ fci 



2k x 



1 

2^ 



(38) 

(39) 

(40) 
fj(l 



where k\ = n (l + av)(l + v)uo p /2c and H2 
a'v){\ — v)u) p /2c. The factor a' is defined analogous to 
Eq. ([2lT) . replacing n Q by 77. a, /?, 7 and 77 are calculated 
with n and n e taken at lu p /2. 

The amplitudes $ oe and $ eo have expressions analo- 
gous to Eq. {!]), with A° e and A^° given by 



A': 



1 



n Q K 



92y) 2 



A^° = Kn eo + (ocp - a)q lx + a p q 2x 



(41a) 



1 \b 



n Q K 



-Aqix + qixf + 7;(qiy + ?2 y ) 2 



-bqfx - gcfiy + {<& ~ Mix - 9<liy) v > (41b) 



where 



V 



2 

n + fj 



Vp + v 



Tin 



V 



2 ' 
n a - fj 



(42a) 
(42b) 



b = (3 p n /ri p , g = 7 p n /r7 p , b = f3n /fj, g = JUo/fj, and 
9? = |qi 1 2 - In Eqs. (|4*Tj) , terms in a and a' were neglected 
with respect to 1, and v 2 was neglected with respect to 
v. Note that /_t oe and /_t eo have a linear dependence on v, 
differently from /i QO , whose dependence is on v 2 . 
In terms of the output angles ^ and £ 2 , we have 



NG{v)£ 



K Kv 

T (Ci + € 3 ) + -2-(€i-€ a ) 



x smc 



/oe(€i,€ 2 » I/ ) 



X exp [-iif« c / oe (^i,^2> I/ )] > 

where JV is a normalization constant and 



(43) 



/oe = Moe + + f)£ 



a 



l.r 



(1 - v)&x 

2 1 = t2 



■f[(l + ^lv + (1-^)6/}- (44) 



$ eo has an expression similar to (|43|) . with 



^-{(i + ^)(KL+5^) + (i-^)e 



4n G 



-5[(l + ^lx + (l-^)6x] 2 

-f [(l + ^fiv + Cl-")^] 2 }- 



(45) 



Fig. [T2] shows density plots of single counts obtained 
from Eqs. (|43H45[) integrated in one of the output angles 
and in v over a bandwidth of 0.025 (corresponding to 
a lOnm filter) for a BBO crystal pumped by a 407nm 
laser, cut for 9 = 42.5°, with L z = 1mm (left) and L z — 
0.25mm (right). Both plots are in good agreement with 
experimental data reported by Lee et al. 
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FIG. 12: Density plot of single counts obtained from Eq. (|47p 
integrated in one of the output angles and in v over a band- 
width of 0.025 (which corresponds to a lOnm filter) for a BBO 
crystal pumped by a 407nm laser, cut for 8 = 42.5°, with 
L z — 1mm (left) and L z — 0.25mm (right). 

The linear dependence of fi oe and fi eo with v leads 
to a frequency spectrum of the two-photon state quite 
different from the type I case. Figure [TBI shows the down- 
converted wavelengths for collinear (^=£2=0) type II 
phase matching as a function of 59, the deviation from 
9 m = 49.223° for BBO with L z = 5mm and L z = 2mm, 
pumped with A p = 351nm. Both plots were obtained 
from Eqs [[331135]) . 



Fig. [TJ] shows the Hong-Ou-Mandel dip for a 1mm- 
long BBO crystal cut for beamlike type II phase match- 
ing [6 — 48.34°), pumped by a 351nm laser. The solid 
line corresponds to the case in which interference filters 
of 20nm bandwidth are used, in good agreement with the 
data reported by Kim [37j • The dashed line corresponds 
to the case in which no interference filters are used. No- 
tice the difference between Figs. [7] and [T3] In type II 
phase matching, the Hong-Ou-Mandel dip is much larger 
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detection Pd is proportional to f dv(\<& oe \ 2 d + |$ eo | 2 ), 
where 



I* 



oe\d 



1$, 



\NG(v)£(0)\ 2 
-KL 



\NG(v)£(0)\ 2 

„:„„2 



KL, , 



oe - ltd - C oe\ 2 )\ , (47a) 
U-CeoH, (47b) 



FIG. 13: (Color online) Down-converted wavelengths for 
collinear (£- L — f 2 = 0) type II phase matching as a function of 
59, the deviation from 9 m = 49.223° for BBO with L z = 5mm 
(left) and L z = 2mm (right), pumped with \ p = 351nm. The 
solid lines correspond to the conditions /i oe =0 and fj, eo =0. 



R n 



R„„ — 



2n fJ, e + n 2 a 2 



2n (i eo + n 2 a 2 
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FIG. 14: (Color online) Hong-Ou-Mandel dip for a lmm-long 
BBO crystal cut for type II phase matching, with 6 = 48.34° 
pumped by a 351nm laser, calculated with Eqs. (|28p and 
(|43|) . The solid line corresponds to the case in which inter- 
ference filters of 20nm bandwidth are used. The dashed line 
corresponds to the case in which no interference filters are 
used. 



than in type I, indicating a narrower frequency spectrum 
in type II. 

It is also possible to write $ oe and $ eo in terms of 
the variables £ s and £ d defined in Eq. [30l but the exact 
expressions are too long to be of practical use. However, 
when 6, 6, g and g are all approximated by 1, we arrive 
at the useful expressions 



F, 



oe — Moe 



*P 7T- a ) 6> 



a P — " ) is 

l + v 



1 c £d 
I I Zn 



(46a) 



(46b) 



When the detectors are scanned in opposite directions in 
the far field, so that £ s = 0, the probability of coincidence 



n a a 



2 + " 



2 + V 



( 2e 




V a 





fOp_ 









x, 



(49a) 
(49b) 



and x is the unit vector in the x direction. 

Expressions (147|) describe two circular profiles, cen- 
tered at £ d = c oe and £ d = c eo , with radii R oe and R eo , 
respectively. Collinear down-conversion occurs when 
Roe = Coe and R eo = c eo . The so-called beamlike type 
II down-conversion occurs in the degenerate case (v = 0) 
when 8 is such that R oe = R eo = 0. 

Let us analyze the case when the detectors are scanned 
in the same direction in the far field, so that £ d = 0. 
Then, 



1$, 



\NG{v)E{KO? 
. 2 f KL * r 

x sine < 



I 2 

\NG{v)£{K£ 
i r KL Z 



x sine 



l^oe 
)i 2 



+ («i.-|)u]}i(50a) 
}.(50b) 



Note that in type II the transfer of angular spectrum 
described by Eqs. ([50)1 is affected by the detuning v 
through the parameters fi oe and (i eo . The detuning has 
the effect of displacing the sine function by an amount 
(n — fj)KL z u/4 in the x axis, in opposite directions for 
|<E> oe |s and I'&eols- In the limit of long crystals, the sine 2 
functions are narrow enough to allow us to write the co- 
incidence detection probability as 



du(\^ oe \ 2 + \^ eo \ 2 s ) 



\N£(K£ S 



\ to / V to / 



,(51) 



where to = (2a p — a)/(n — rj). In type II phase match- 
ing, due to the linear dependence of y, oe and y, eo on v, 
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FIG. 15: (Color online) Comparison between experimental re- 
sults (•) and the prediction of Eq. (|50[) (solid line) for a BBO 
crystal with L z = 5mm in collinear type II phase match- 
ing, pumped by a 405nm laser with wq ~ 25/im. The two 
detectors are scanned in the same sense, along the x direc- 
tion. The Gaussian predicted by Eq. (|50[) coincides with the 
expected profile of the pump beam. The dashed line shows 
the expected curve for the monochromatic case. 

the angular spectrum of the pump field may be trans- 
ferred to the two-photon state if the bandwidth of the 
detection filters is broad enough, even for long crystals. 
Figs. [T5l and [TBI show comparison of the predictions of Eq. 
(I50p with experimental data for a 5mm-long BBO crystal 
pumped by a 405nm laser with wo ~ 25/im and detection 
filters with a bandwidth of lOnm. In both figures, there 
is a good agreement between experimental data and the 
expected Gaussian profile of the pump beam. The ex- 
perimental conditions are the same of Figs. [10] and [11] 
except for the phase matching type. Notice the difference 
between the two cases when the detectors are scanned 
in the x direction (Figs. [10] and [T5]) . In the case of 
Fig. [T5] the lOnm filter allows a complete transfer of 
the pump beam profile to the coincidence detection. The 
dashed line shows the expected curve for the monochro- 
matic case, that is, G{v) as 8{v) (zero bandwidth filter). 



VI. SUMMARY AND CONCLUSION 

We have discussed in detail the two-photon state gen- 
erated by spontaneous parametric down-conversion in 
bulk crystals, taking into account the effect of crystal 
anisotropy, in both in type I and type II phase match- 
ing. Our discussion was based on the perturbative ap- 



proach introduced by Hong and Mandel, in the context of 
Fourier Optics, in which the fields are treated by means of 
their plane-wave expansion. The frequency spectrum of 
the down-converted fields was also considered, in the ap- 
proximation of small detuning and monochromatic pump 
beams. Extension of the theory to cover pulsed pump 
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FIG. 16: (Color online) Comparison between experimental re- 
sults (•) and the prediction of Eq. (|50|l (solid line) for a BBO 
crystal with L z = 5mm in collinear type II phase match- 
ing, pumped by a 405nm laser with wo ~ 25/^m. The two 
detectors are scanned in the same sense, along the y direc- 
tion. The Gaussian predicted by Eq. (|50[) coincides with the 
expected profile of the pump beam. 



beams seems straightforward. Several approximations 
were made in order to provide simple expressions that 
still exhibit the main features of the frequency and spa- 
tial spectral properties of the two-photon states. These 
expressions were shown to be in good agreement with 
experimental data for some selected situations. In par- 
ticular, the spectral content of the two-photon state and 
the transfer of angular spectrum from the pump beam 
to the two-photon state were analyzed, and a significant 
differences of these features in type I and type II phase 
matching were discussed. The results presented here are 
a contribution to the Fourier optics of two-photon states 
and may be helpful to improve the understanding and 
further development of entangled state sources for quan- 
tum information and quantum communication. 
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